Manipulating Majorana fermions using supercurrents 
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Topological insulator edges and spin-orbit-coupled quantum wires in proximity to s-wave super- 
conductors can be tuned through a topological quantum phase transition by a Zeeman field. Here 
we show that a supercurrent flowing in the s-wave superconductor also drives such a transition. 
We propose to use this mechanism to generate and manipulate Majorana fermions that localize at 
domain walls between topological and nontopological regions of an edge or wire. In quantum wires, 
this method carries the added benefit that a supercurrent reduces the critical Zeeman field at which 
the topological phase appears. 
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Introduction. — Emergent Majorana fermions in a con- 
densed matter setting are currently attracting much 
attention.'Il'^ Zero-energy Majorana fermions comprise 
the simplest non-Abelian anyon and promise fascinat- 
ing applications to topological quantum information 
processing.!^ Recently, the set of candidate systems sup- 
porting Majorana fermions has greatly expanded beyond 
quantum Hall systems^ ^ with the realization that several 
materials can be driven into a topological superconduct- 
ing phase when placed in proximity to a conventional 5- 
wave superconductor. This was initially understood for 
topological insulators followed by 2D 5-wave super- 
fluids with Rashba spi n-orbi t interaction, 1^ spin- orbit- 
coupled quantum wellJ^^I^ and nanowires,'^^'^ half- 
met als,'^^^^ and 3D topological insulator nanoribbons.!^ 

Nanowire proposals are attractive as they involve 
widely available materials and provide detailed recipes 
for manipulating the Majorana fermions^^ — a prerequi- 
site for verifying their non-Abelian statistics and per- 
forming topological quantum information processing. 
While initial proposals examined simple mean-field mod- 
els of clean wires proximate to a superconductor, more 
recent work indicates that the induced topological phase 
persists in the presence of moderate interactions^^^ or 
disorder, ^^ ^^ and c onsid ered setups for probing the Ma- 
jorana bound states .H^^ Experimental challenges nev- 
ertheless remain: realizing the topological phase requires 
control over the wire's global electron density and the 
application of significant Zeeman fields without destroy- 
ing superconductivity. Furthermore, manipulating Majo- 
rana fermions by loc ally con trolling the electron density 
using gate electro des^^l^II^ is nontrivial due to strong 
screening by the superconductor. 

Here we show that the latter two challenges can be 
greatly alleviated by applying supercurrents in the bulk 
superconductor. These supercurrents cause a spatial gra- 
dient of the phase of the proximity-induced pair poten- 
tial in the wire, which drives a transition between the 
nontopological and topological superconducting phases. 
Remarkably, the supercurrent also allows one to access 
the topological phase at weaker Zeeman fields. Spatially 
varying the phase gradient along the wire moreover gen- 



erates controllable Major ana-carrying domain walls be- 
tween nontopological and topological regions. Switching 
the supercurrents (and hence the phase gradient) along 
sections of the quantum wire transports these Majorana 
fermions preserving the gap, obviating the need for local 
gating. Schematically, this can be achieved by the device 
shown in Fig. [l] Our scheme, in fact, applies equally 
well to the edge of a two-dimensional topologi cal in su- 
lator with proximity-induced superconductivity.H'^ We 
start by elucidating the physics in this setting, since it 
is somewhat simpler to analyze, and then turn to the 
quantum wire case. 

Topological insulators. — Our analysis begins from the 
Bogoliubov-de Gennes equation describing a topological 
insulator edge proximate to an g-wav e superconductor 
and subjected to a magnetic 

fieldPini 

n = (^/pcr^-/i)r^-5a^ + Ae^^^^V+ + Ae-^^(^V_. (1) 

Here, p is the momentum along the edge, u measures 
the edge-state velocity, ji is the chemical potential, A 
and (j){x) denote the magnitude and (position-dependent) 
phase of the proximity-induced pair potential, and 5 > 
is the Zeeman field. The Pauli matrices and respec- 
tively act in spin and particle- hole space. Equation ([T]) 
is written in a Nambu basis with spinors of the form 

In the absence of a phase gradient, say (/) = 0, the 
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FIG. 1. Schematic device for manipulation of Majo- 
rana fermions by supercurrents. Different sections of the 
proximity-inducing s-wave superconductor carry different 
currents, putting the quantum wire into an alternation of 
topological (T) and nontopological (NT) phases. Changing 
the supercurrent in the segments allows for manipulation of 
the Majorana fermions 7 associated with the domain walls. 
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FIG. 2. Phase diagram for the topological insulator case with 
/i = 0. Increasing the Zeeman field B or the superconducting 
phase gradient V(/) induces a phase transition from the topo- 
logical phase (TP) into a nontopological phase (NTP). The 
latter has two regimes, separated by the dotted line, depend- 
ing on whether the minimal excitation energy Amin is located 
at zero or nonzero p. The red dots indicate representative 
points at which the spectrum is illustrated in Fig. [3] 

spectrum of Eq. ([T]) is easily derived by repeatedly squar- 
ing exploiting the fact that the Bogoliubov spectrum 
is symmetric about zero energyP^ Focusing for simplic- 
ity on /i = until specified otherwise, one finds energies 
E±{p) = {{up)^^\A±B\^y/^. With B < A,n describes 
a topological superconducting state^^ similar to that of 
Kitaev's model for a ID p-wave superconductor.^^ The 
gap closes when 5 = A, signifying a topological quantum 
phase transition into a trivial superconducting state, and 
reopens at larger B; see Fig. [2j 

Consider now Eq. ([T]) with a nonuniform phase (j){x). 
Rather than studying the Hamiltonian ([T]) directly, it is 
advantageous to gauge away <p{x) from the pairing term 
using the unitary transformation 



U = exp{i(l){x)Tz/2}. 



(2) 



Note that this gauge transformation multiplies the elec- 
tron and hole components of the Bogoliubov-de Gennes 
spinor by opposite phase factors. The transformed 
Hamiltonian W = UHU^ becomes 



H' 



u I p 



V(t) 



Tz-B(Jz^At^, (3) 



which depends only on the gradient of ^, making gauge 
invariance explicit. 

For uniform phase gradients along the wire, the 
Hamiltonian %' in Eq. ([3| is readily diagonalized; Fig. |2] 
summarizes the resulting phase diagram. Representative 
spectra, corresponding to the points in the phase diagram 
marked in Fig. [2j appear in Fig. |3] When the system 
resides in the nontopological phase {B > A) one finds an 
enhancement of the gap as the phase gradient increases 
from zero. Conversely, when the system begins in the 
topological phase {B < A), the gap initially decreases 
when applying a nonzero phase gradient. Eventually the 
gap closes at the critical phase gradient 



(V</.)* = (2A/«)[l-(i?/A)^] 



211/2 



(4) 



which signals the transition into a nontopological phase 
as indicated by the full line in Fig. [2j Note that (V(/))* 
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FIG. 3. Bogoliubov de Gennes spectrum of the topological 
insulator at /i = for various Zeeman fields and supercon- 
ducting phase gradients, as indicated in the figure. The cor- 
responding points are also marked in the phase diagram of 
Fig. |2] By symmetry, only the positive-energy dispersion is 
shown. E and p are measured in units of eso and 2mu, re- 
spectively. 



is of the order of the inverse of the proximity-induced 
(zeio-B) coherence length of the edge state, ^ = u/A. 
Increasing the phase gradient beyond (V^)* reopens the 
gap for any nonzero B. (At zero 5, there is a transition 
into a gapless phase.) This shows that a gradient of the 
phase of the gap function, originating from supercurrents 
in the bulk 5- wave superconductor, can indeed induce 
a transition between the topological and nontopological 
superconducting phase. One can understand this result 
intuitively by observing from Eq. (|3| that the phase gra- 
dient behaves exactly as a magnetic field oriented along 
X, which effectively shifts the electrons' momentum (in 
contrast to the Zeeman field B which opens a gap). When 
= and A > 5, electrons with opposite momenta 
are resonant and can easily form Cooper pairs, driving 
the system to the topological phase with a A-dominated 
gap. A nonzero phase gradient, however, breaks the res- 
onance between p and —p states and suppresses Cooper 
pairing. If this shift is too large, the Zeeman field B 
dominates the gap and the edge forms a trivial phase. 

When increasing the phase gradient even further in the 
nontopological phase, there is another characteristic line 
where the minimal gap in the excitation spectrum is no 
longer located at p = 0, but rather at a finite p. This 
crossover line is indicated by the dotted line in Fig. [2] 

Quantum wire. — We now turn to the phase diagram 
for a quantum wire, proximity coupled to a supercurrent- 
carrying s-wave superconductor. Just as for the topolog- 
ical insulator edge, this phase diagr am is obtained from 
the quantum wire HamiltoniaiP^EH 



n = (pV2m+^/pcr^-/i)r^-5cr^+Ae^^(^V++Ae" 



-i(f){x) 



(5) 



by performing the gauge transformation In this 

Hamiltonian, m denotes the effective mass and u the 
Rashba spin-orbit-coupling strength. When (j){x) is uni- 
form, the quantum wire forms a nontopological supercon- 
ducting phase for 5 < A, undergoes a topological quan- 
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FIG. 4. Phase diagram of a quantum wire proximity coupled 
to a supercurrent-carrying s-wave superconductor for A = Cso 
in the Zeeman field-phase gradient plane at /x = (main 
panel) and the chemical potential-phase gradient plane at 
5/ A = 1.2 (inset). In both panels the full (dashed) line 
indicates where the gap vanishes at zero momentum (near 
p = ±Pf), leading to gapped topological (TP) and nontopo- 
logical (NTP) phases as well as a gapless phase (GLP). The 
red dots indicate representative points at which the spectrum 
is illustrated in Fig. [5] 
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FIG. 5. Bogoliubov de Gennes spectrumof the quantum 
wire in proximity to a supercurrent-carrying s-wave super- 
conductor, focusing on positive energies due to the symmetry 
(PjE) (— p, —E). The chemical potential has been set to 
/i = 0; othe parameters are specified above and indicated by 
red dots in Fig. |4] E and p are measured in units of eso and 
2mu, respectively. 



turn phase transition at 5 = A, and enters a topological 
superconducting state at larger Zeeman fields. 

To understand the resulting phases when a uniform 
phase gradient is present, we first focus on the excitation 
spectra for /i = in the vicinity of the Fermi points (for 
A = 5 = 0)atp = and ±pf ±2mu. Figure [I] 
shows the resulting phase diagram in the Zeeman-field- 
phase gradient plane. The full line indicates where the 
gap in the excitation spectrum closes at p = while re- 
maining finite on both sides of the line. The dashed line 
meanwhile indicates where the gap closes near the Fermi 
momentum pp. In this case, the gap is finite only on the 
low phase gradient side of the line, while it remains closed 
on the high phase gradient side. Figure |5] shows repre- 
sentative excitation spectra for the parameters indicated 
by red dots in the phase diagram of Fig.ga). These fig- 
ures show that for B below a critical value 5* a uniform 
phase gradient drives the wire from the nontopological 
phase directly into a gapless state. More interestingly, at 
5* < 5 < A we find that the quantum wire enters the 
topological phase at intermediate phase gradients before 
reaching the gapless state at even larger phase gradients. 
Finally, when B > A the quantum wire forms a topo- 
logical phase all the way up to a critical phase gradient 
where it becomes gapless. 

Quantitatively, we find for /i = that the phase bound- 
ary between the gapped topological and nontopological 
phases (solid line in Fig. [4]) is implicitly given by 



4eso 



^i(V0)* 



4e: 



so 



B^ 



(6) 



where ego = mv? /2. Note that for \B'^ — /S?\ 
this yields the phase boundary ix(V^)* = 2{/S? - B'^Y^'^ 
which is independent of ego- For A <C Cgo, the junction 
of the full and the dashed lines in Fig. |4] can be readily 
accessed analytically and occurs at 5* = A^/4eso. 



The inset in Fig. |4] illustrates the effect of a phase gra- 
dient at a nonzero chemical potential for 5 > A. At 
= the wire is topological for < — A^. 

Remarkably, a finite (but not too large) phase gradient 
extends the topological region to larger making the 
existence of Majorana bound states less sensitive to the 
tuning of the chemical potential. 

As for the topological insulator edge, these phase dia- 
grams indicate that a phase gradient induces a topologi- 
cal phase transition over a wide region of Zeeman fields, 
allowing one to realize Major ana-carrying domain walls 
by spatially varying supercurrents. Moreover, these re- 
sults show that in the presence of a phase gradient, the 
Majorana carrying topological phase can be induced in 
the quantum wire at weaker Zeeman fields. Both of these 
facts may significantly simplify the experimental realiza- 
tion of Majorana fermions in quantum wires. 

To further illustrate the supercurrent-driven transition 
into a topological phase, we computed the lowest two ex- 
citation energies of a finite-length wire subjected to a 
uniform phase gradient across the entire system using 
a lattice model that recovers Eq. ([5| in the low-density 
limit. Fi gure [g] shows the results for a 1000-site wire 
with /i = 0, B/A = 0.9, and A/ego = 2. Since B < A 
here, one finds that both excitation energies remain fi- 
nite up to a critical phase gradient, indicating that the 
system forms a nontopological phase. Above the critical 
phase gradient, the lowest excitation energy drops essen- 
tially to zero due the formation of localized Majorana 
end-states associated with the entry into the topological 
phase, while the second energy remains finite reflecting 
the wire's bulk gap. Beyond a second critical phase gra- 
dient, the bulk gap closes and the system enters a gap- 
less phase, in agreement with the phase diagram in Fig. 
|4] The lower left of Fig. [6] displays the probability dis- 
tribution for the near-zero mode generated by a phase 
gradient uV(j) /{2A) = 0.6 induced only over the central 
half of the same 1000-site wire. The phase gradient ere- 
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FIG. 6. Lowest two excitation energies of a 1000-site quantum 
wire vs. the phase twist generated by a supercurrent through 
the entire wire for /x = 0, B/A = 0.9, A/eso = 2. At zero 
twist, the wire is in the nontopological phase with B < A. 
At larger twists, the wire enters the gapped topological phase 
where the lowest excitation energy drops to zero — due to the 
formation of Majorana modes — before reaching the gap less 
phase at even higher phase gradients, in line with the phase di- 
agram in Fig. [4] Inset (lower left): Probability distribution for 
the near-zero mode obtained when current flows only through 
the central half of the wire, creating a trivial-topological- 
trivial domain structure with two Majorana modes. 



ates two domain walls between the trivial outer ends and 
topological inner region, each of which clearly binds a 



localized Majorana mode as expected. 

Estimates. — The proposed scheme to manipulate Ma- 
jorana fermions requires one to establish a sufficient 
phase gradient to drive the quantum wire (or topological 
insulator edge) between the topological and nontopolog- 
ical phase without hitting the critical current density of 
the s-wave superconductor. This critical current den- 
sity corresponds to approximately one 27r phase wind- 
ing within the superconductor's coherence length ^sc- In 
comparison, the topological phase transition occurs for a 
phase gradient of roughly one phase winding per coher- 
ence length of the wire ^wire = u/A. Thus, our scheme 
requires ^ ^wire-^^ Both parameters entering the co- 
herence length ^sc = ^f/Asc of a clean superconductor, 
namely the Fermi velocity vp and the pairing amplitude 
Age, typically exceed the corresponding quantities for the 
wire. Thus, if both the wire and superconductor are 
clean, it is in general possible to satisfy the requirement 
by a suitable choice of materials and device parameters. 
The estimates become yet more favorable for a dirty 5- 
wave superconductor with coherence length ^ = (^sc^)^^^, 
where £ is the elastic mean free path. 
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of the smaller required Zeeman field is not compensated by 
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nitude A of the proximity induced gap due to the applied 
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